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Abstract 

Despite the abundance and importance of three-dimensional systems, relatively little progress 
has been made on spatiotemporal nonlinear optical waves compared to time-only or space-only sys¬ 
tems. Here we study radiation emitted by three-dimensionally evolving nonlinear optical waves in 
multimode fiber. Spatiotemporal oscillations of solitons in the fiber generate multimode dispersive 
wave sidebands over an ultrabroadband spectral range. This work suggests routes to multipurpose 
sources of coherent electromagnetic waves, with unprecedented wavelength coverage. 
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Eigenmodes are ubiquitous tools for describing complex wave systems. For nonlinear com¬ 
plex wave systems, however, the superposition principle is not applicable. In special cases, 
nonlinear wave systems possess solitons, which act to some extent as nonlinear eigenmodes. 
Combined with more general nonlinear attractors and perhaps insights from linearized sys¬ 
tems, researchers may thus build up a conceptual understanding of complex nonlinear wave 
dynamics. In optics, one-dimensional (ID) dynamics in single-mode waveguides have been 
thoroughly explored, with many advances hinging on the robust nonlinear attraction of 
solitons [IH3]. In unbounded 3D systems, dynamics have been successfully explained largely 
in terms of a nonlinearly attracting instability, namely spatial or spatiotemporal collapse |1]. 
In reality, the collapse singularity is avoided by a multitude of higher-order effects, and 
the field eventually expands, by propagating linearly. Post-collapse, promising results have 
been obtained considering conical waves, which are the eigenmodes of the 3D linear wave 
equation^ |6]. Although not yet widely adopted, conical wave solutions to the nonlin¬ 
ear wave equation may provide even deeper insight [71 E]. Given the great advantages the 
concepts of solitons and collapse have provided in the field of nonlinear optical waves for 
studying single-mode waveguides and free space filamentation, it is natural to seek simi¬ 
larly useful nonlinear objects in multimode waveguides. Multimode waveguides include as 
limiting cases single mode (ID) and free-space (3D), so these hypothetical objects, whether 
solitons, nonlinearly attracting instabilities, guided conical waves, or something else entirely, 
could even help to conceptually unify nonlinear optical dynamics across dimensions. 

Solitons in optical single-mode fiber (SMF) have been intensely researched, because they 
are relatively accessible both analytically and experimentally. Equally important, soliton 
dynamics are critical to telecommunications, mode-locked hber lasers, and compact white- 
light sources with high spatial mode quality. Multimode fibers (MMFs) could provide major 
benehts for various applications, from spatial division multiplexing in commnnications[9l ITO] . 
to high-power, versatile fiber lasers and white-light sources [TT]. Although wave propagation 
in MMF is still experimentally and theoretically challenging, recent theoretical advances |12h 
[17] combined with opportunities for major performance increases provide ample motivation 
for their study. From a scientific perspective, MMF is an ideal environment for studying 
spatiotemporal nonlinear optical dynamics. By judicious design, or by control of the initial 
excitation, researchers may control the spatiotemporal characteristics of nonlinear dynam¬ 
ics, through variation of the effective dimensionality, the coupling between modes, or their 
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individual dispersions. 

These factors have motivated recent work on nonlinear optical waves in MMFs[Tn [T5l 
flGl [T81 - I22] . In particular, control of the spatial excitation of a MMF provides a means 
of controlling the nonlinear spatiotemporal dynamics, and therefore the characteristics of 
a supercontinuum[TT]. Launching «i200-nJ and 500-fs pulses at 1550 nm into a graded- 
index (GRIN) MMF produces remarkable visible light emission characterized by a series of 
spectral peaks with non-uniform spacing. Numerical simulations confirm these peaks, and 
suggest that even more spectacular emission occurs at long wavelengths, where fused sil¬ 
ica is effectively opaque (Figure 1). These observations are puzzling and fascinating: some 
novel nonlinear optical mechanism, unidentified by prior fiber-based or bulk supercontin¬ 
uum research, generates and links electromagnetic radiation over 2 orders of magnitude in 
wavelength (from > 50 to < 0.5 pm). 


Wavelength (nm) 



Frequency (THz) 


FIG. 1. Simulated (top) and example experimental (bottom) supercontinunm in multimode GRIN 
fiber. The pump pulse at 1550 nm creates a spectrum with a series of red-shifted and blue-shifted 
peaks. In the bottom panel, the y-axis reference (0 dB) is the maximum intensity of the 1550 nm 
pump peak. 


Here we provide a theoretical explanation for this phenomenon: dispersive waves are gen¬ 
erated by the spatiotemporal oscillation of multimode solitons. The process is inherently 3D, 
with spatiotemporally-evolving nonlinear waves emitting spatiotemporally-evolving disper- 
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sive waves. However, we show that insights from solitons of the ID nonlinear Schrodinger 
eqnation (NSE), can prove usefnl in nnderstanding the higher-dimensional system. This 
nnderstanding also snggests rontes to generating nltrashort pulses in hber at wavelengths 
outside current capabilities, and may provide a means of interfacing normally distinct parts 
of the electromagnetic spectrum. 

Solitons are solutions to a conservative equation. To use them in applications such as 
mode-locked lasers and optical telecommunications [23] - [3Tj . loss must be compensated. This 
is naturally accomplished by laser gain, or by amplihers placed throughout the line. The 
periodic perturbations of gain and loss, however, can destabilize a soliton. The origin of the 
instability, and the characteristic spectral sidebands that are its signature, is the fact that a 
periodic perturbation can phase-match dispersive wave emission at particular frequencies [351 
[36] . For a perturbation period (spacing of amplihers, or laser cavity length) Zc, the phase¬ 
matching condition is approximately: 

{ksoi - kdis) = 2rmilZc ( 1 ) 


where m is an integer, and kdis is the wavevector of the dispersive wave, kgoi is the soli- 
ton wavevector, which in the absence of higher-order dispersion is equal to |/ 92 |/ 2 r^ = 
7 r/( 4 Zo),where (32 is the group velocity dispersion, r is the soliton duration and Zg = 
7i/{4ksoi) is the soliton period. This quasi-phase matching leads to resonant emission at 
frequencies separated from the pump by 


k^res 


8Z„m 


- 1 


( 2 ) 


. This result was later rehned to include the third-order dispersion of the hber[3ST 
, where resonant frequencies were found to be roots of the equation: 


1 

{ksoi - kdis) = 217171 (Zc = (- + 463 - 863 ) - (—— 




(3) 


where 63 = —f3^/{Q\l32\T) {(3^ is the third-order dispersion) and D is the angular frequency 
separation from the pump (in units of t~^). Production of dispersive waves is the primary 
limitation to the performance of soliton hber lasers. 

Although they were predicted as early as the 198ns [TT1 - l45] . solitons consisting of pulses 
within multiple spatial modes have only recently been studied experimentally [TTl [T 8 l [ 22 ]. 
Initial work shows that solitary waves, termed MM solitons, can be excited in multimode 
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graded-index (GRIN) fibers. Similarly to single-mode solitons, these pulses resist group- 
velocity dispersion and adjust themselves adiabatically in response to perturbations such as 
intrapulse Raman scattering. This behavior makes them a useful conceptual tool for under¬ 
standing complex nonlinear processes beyond the NSE description, such as those involved in 
supercontinuum generation [TTl 1^ . Unlike single-mode solitons however, MM solitons addi¬ 
tionally resist modal velocity dispersion - he., the fact that each of the numerous consitutent 
modes of the soliton has a different group velocity. Fission of MM solitons is spatiotem- 
poral: it yields multiple MM solitons which can have many different modal distributions. 
In fact, MM solitons can have a bewildering range of spatiotemporal shapes, which makes 
their dynamics much richer and challenging than single-mode solitons. While they are in 
some sense natural extensions of the ID NSE soliton to the case of a MM fiber, they prob¬ 
ably do not fulfill the most rigorous definitions of 'soliton' (some special cases are already 
known, however [T^ [16] ). Questions such as whether they are stable over very long distances, 
how they interact with one another, and how many modes can be involved, remain largely 
unanswered. 

When a beam excites multiple modes of a GRIN fiber, it propagates through the fiber 
with a characteristic spatial oscillation with period (pitch) P = ttR/ a/ 2A, where R is the 
core radius and A is the core-cladding index difference [IHl SO]- This oscillation causes the 
intensity of the beam to periodically evolve, as in a loss-managed soliton transmission line 
or soliton fiber laser. For a pulsed beam, oscillations will only occur as long as the pulses in 
each mode maintain co-localization and a phase relationship [18]. Hence, a soliton containing 
multiple spatial modes experiences a periodic oscillation of its peak intensity and therefore 
is likely to emit dispersive radiation at particular frequencies depending on the period of the 
oscillation. 

We consider a simulation using the generalized multimode nonlinear Schrodinger equa¬ 
tion (GMMNLSE)[T2|, with the first five radially-symmetric modes excited uniformly for 
simplicity (Figure 1). Experiments are conducted with the procedures described in Ref. 
HD- In both simulation and experiment, we find that the soliton oscillation-induced dis¬ 
persive waves (ODWs) are observed for many different initial spatial conditions, including 
non-radially-symmetric excitations. Typically we observe that the ODWs are more appar¬ 
ent when the intensity of the initial pulse is higher (either by large energy, tight spatial 
localization, or both). Ghanging the initial noise level or the pulse duration results in dif- 
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FIG. 2. Dynamics of dispersive wave formation in GRIN fiber, a) Temporal and spatial breathing 
of the field (inset: zoom in near the onset of dispersive wave generation), b) Evolution of the 
spectral intensity of the whole field through the same distance, c) Energy in each dispersive wave 
band. The x-axis scales are normalized to the linear spatial oscillation period of the GRIN fiber, 
equal to 407 fim. These dynamics are also shown in Supplementary Movies I and 2, which provide 
a considerably more complete representation of the complex spatiotemporal evolution. 


ferent dynamics. Nonetheless, provided the pulse energy and fiber length are sufficient, and 
that multiple modes are excited, the qualitative features (including the spectral positions) 
of the ODW emission process are similar. In simulations, we observe that the energy of each 
ODW is distributed roughly equally among all the modes, with the red-shifted (blue-shifted) 
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FIG. 3. a. ODW from MMGNLSE simulations (MMGNLSE), eqn. 4 with 7 ( 2 ;) as in text,ID NEE 
simulations, and an experimental example spectrum, b. Eunctions used in the quasi-lD approxima¬ 
tion, along with the peak intensity variation of the field in MMGNLSE simnlation. c.Comparison 
of the periodic ID phase-matching model with multimode simulations by the GMMNLSE and ex¬ 
periments in GRIN MMEs. Continuous curves are values of m in Eq. 3, plotted with the best-fit 
parameter values. 

sidebands exhibiting a slight preference for the low-order (higher-order) modes. 

Representative simulation results for a 400-fs, 1550-nm MM pulse with energy 164 nJ 
launched into a GRIN MMF are shown in Figure 2. (The results are also shown as Sup¬ 
plementary Movies 1 and 2.) The spatial and temporal breathing of the held are evident 
in Fig. 2a. While we show only the blue-shifted ODWs (Figure 2b), the simulation also 
predicts red-shifted ODWs. These are remarkably outside the transparency window of fused 
silica (the hrst appears at roughly 72 THz (4200 nm)). As the pulse traverses the hber, 


7 























the spectrum develops a dispersive wave near 300 THz (1000 nm), as well as ODWs in the 
visible and mid-IR regimes. Attenuation is included in the simulations with an assumed fre¬ 
quency dependence a = 01550 exp—(/ — c/(1550nm))/6;, where 01550 = O.OSdS/fcm is the 
attenuation (units m“^) at 1550 nm, and bi = 0.0062 PHz 80 nm) models the increasing 
loss into the infrared. For these parameters, attenuation is ~ 0.1 MdB/km at 4200 nm. We 
neglect loss above the 1550-nm base level on the blue side of the spectrum. Despite this 
enormous loss, appreciable energy is still produced in the mid-IR ODWs. This is partly 
due to the relatively short length of the hber (150 periods correspond to only 6 cm), and 
to the enormous gain the ODWs experience (relative to the initial continuum level). The 
energy in each ODW may grow superexponentially, with the hrst blueshifted (redshifted) 
ODW experiencing roughly 70 dB (50 dB) net gain in 1 cm (Figure 2c). After about 100 
oscillation periods the ODW energy in each spectral band saturates to a nearly constant 
value. The observed superexponential growth and subsequent saturation of the ODW energy 
is in qualitative agreement with the analytical prediction of the soliton perturbation theory 
that was developed to describe large-amplitude soliton intensity oscillations in periodically 
amplihed transmission hnks[171 HE]. On the other hand, for small amplitude oscillations of 
the soliton, the theory only predicts a linear growth of the ODW with distance |49]. 

To verify that the oscillations underly the ODWs, we add oscillations artihcially to a 
single-mode NSE. This is done by making the nonlinear coefficient a periodic function of the 
longitudinal coordinate, so that the pulse evolution equation is 


aAu,t) ^ _A d^A(zj) ^ 


(4) 


dz 2 dt"^ 6 dt'^ 

where A{z,t) is the pulse envelope, and 'y{z) is the ; 2 -dependent nonlinear coefficient. Fig¬ 
ure 3 compares the result of solving this equation to the results above found using the 
GMMNLSE, and with experiment. Figure 3a shows the solution of equation 4 with the 
indicated form of 'y{z). Due to the sinusoidal oscillation of the beam radius, one expects an 
accurate 7(2;) to be of the form 7(2;) = ^o/(,[rmax - Tmean] + Tmean sin {2'kz/P)Y. We use the 
RMS widths for rmax and r^ean here, and rescale 70 appropriately, since the usual spatial 
width measurement (mode held diameter) is not well-dehned for the complex spatiotemporal 
helds that occur. At the onset of the ODW generation («i 145-155 periods), Vmax = 14.5 /rm 
and Tmin = 5.2 pm, so that 7(2:) = 7o/(9.8-|-4.6sin (27rz/P))^ (IDb, Fig. 3c). The observed 
intensity oscillations of the MMF held (MMGNLSE, Fig. 3b) are approximated better by 





7o/(9.8 + 6.9 sin (27r2:/P))^ (space+time fit, Fig. 3b). This is because the MM soliton’s 
duration also oscillates: the soliton tends to lengthen when its spatial width increases and 
shorten when its spatial width decreases (Fig. 2a and inset, Supplementary Movie 1). It is 
this spatiotemporal oscillation that generates the ODWs. 

The experimental peak locations are consistent with simulation and analytic theory. Fig¬ 
ure 3c shows the results of fitting the experimentally-measured and simulated peak locations 
with the roots of Eq. 3 for various m. For the experimental (simulated) peaks, fitting yields 
(32 = -26 (-25) fs^/mm, and /33=143 (143) fs^/mm. In both cases, the peaks are fit by the 
simple theoretical approximation well, with the discrepancy attributed to the slightly differ¬ 
ent dispersion characteristics of the real and simulated hbers. The optimal values are near 
the simulation parameters, except that the optimal (32 is greater by 16%, which suggests 
that modal dispersion makes an effective contribution. 

Although the measured ODW positions are well-described, their amplitudes can vary 
due to several effects beyond both the simulations and analytic approach. First, due to the 
presence of many modes in the hber, phase-matched intramodal four-wave mixing (FWM) 
may occur involving the 1550-nm pump, its primary dispersive wave at 800-1000 nm, various 
ODWs, and Raman-shifted MM solitons. For the low-order ODWs in particular, we note 
that there are several energy-conserving FWM processes. These may amplify red-shifted 
ODWs at the expense of specific blue-shifted ODWs. Particularly in the presence of many 
modes, intramodal FWM can be a complex process |50], particularly when cascaded mixing 
is considered [51]. It is probably why, for certain initial spatial conditions, various low-order 
blueshifted ODWs are attenuated or even invisible experimentally. We obtained the peak 
locations in Figure 3c from experiments with multiple initial conditions, in order to account 
for low-order blue-shifted ODWs that were attenuated for any given intial condition. Sec¬ 
ond, because the spacing of the ODWs is quite close to the Raman bandwidth of fused silica, 
certain ODWs may experience Raman gain from one another and from the third-harmonic 
light (THG: its efficiency may also benefit from intramodal phase-matching). Lastly, the dy¬ 
namic range and spectral resolution of the spectrometer limits the visibility of low-amplitude 
features, and broadens narrow spectral features. 

Another remaining mystery is the relatively high amplitude of the ODWs observed in 
experiments, compared to simulations. Intramodal FWM may play some role, as well as 
larger oscillations and THG. Figure 3a shows the spectrum obtained by solving the ID 
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nonlinear envelope equation (NEE) with j(z) = 7o/(15.5 + 14.5 sin (27rz/P))^, including 
self-steepening and assuming averaging of the oscillating Raman integral[n2]- Larger oscilla¬ 
tions produce relatively more intense dispersive waves, because the soliton is more strongly 
perturbed. Larger oscillations occur when more modes are coherently locked together |22]. 
and as the experiment contains much more than 5 modes, we choose a functional form of 
7 (z) to model this. In addition, THG provides an initial seed in the visible, leading to signif¬ 
icantly enhanced resonant emission at nearby wavelengths. Periodic backconversion occurs 
for frequencies near 3uo, so the largest enhancement is for slightly smaller frequencies. 

The ODWs generated by soliton resonances in MMF have relevance to applications. 
For example, by hltering the hrst redshifted or blueshifted ODW, one may generate pulses 
in wavelength regions well outside the gain spectrum of available hber dopants. Tuning 
may be achieved by changing the pump wavelength or hber pitch. In fact, because the 
modulation instability gain spectrum of a CW held at the pump wavelength overlaps with 
the soliton sidebandspS], hltered ODWs could be parametrically amplihed (using either the 
circulating pump pulse or an injected CW held) provided some mechanism was introduced 
to compensate for the chromatic walk-oh between the pump and ODW helds. ODWs spaced 
over a very broad wavelength range are generated coherently from the optical pump pulse. 
Given the remarkable range and coherence, spatiotemporal soliton oscillation may provide a 
means of generating synchronized ultrashort pulses in diherent regions of the electromagnetic 
spectrum. In an appropriate waveguide, an optically-derived coherent connection between 
microwave, deep ultraviolet and optical frequencies may be possible. 

In summary, we have shown that the spatiotemporal oscillation of nonlinear waves in 
GRIN multimode hber causes the generation of their spatiotemporal dispersive waves. These 
dispersive waves can be described relatively well by simulations using the GMMNLSE, 
and insight can be gained by approximating the dynamics in a quasi-ID model with a 
longitudinally-varying nonlinearity. Future work, involving more advanced models and ex¬ 
perimental methods, can answer a few of the open mysteries about this process, including 
verihcation of these hypotheses, the spatiotemporal dynamics leading to the multimode 
supercontinuum, and the spatiotemporal structure of the dispersive waves. This work illus¬ 
trates valuable conceptual connections between ID and the intermediate dimensional system 
of the GRIN MMF. Practically, it provides a route to hber-based ultrashort pulse sources 
with tunable wavelengths far outside the range of any current hber-optic technique. 
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